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IN-MEDIUM ION MASS RENORMALIZATION AND LATTICE 
VIBRATIONS IN THE NEUTRON STAR CRUST 



The inner crust of a neutron star consists of nuclei immersed in a superfluid neutron 
liquid. As these nuclei move through the fermionic medium they bring it into motion 
as well. As a result their mass is strongly renormalized and the spectrum of the ion 
lattice vibrations is significantly affected. Consequently we predict that the specific heat 
and the lattice thermal energy of the Coulomb crystal at these densities are noticeably 
modified. 

In the standard picture the inner crust of a neutron star consists of neutrons, 
protons and electrons. At low temperatures ( T < 1 MeV) one expects the pres- 
ence of atomic nuclei (ions) immersed in a neutron gas. At the densities: 10 11 
g/cm 3 < p < 10 14 g/cm 3 they form a crystal lattice stabilized by the Coulomb in- 
teraction. The electrons which are ultrarelativistic at these densities form a strongly 
degenerate, uniform gas. The detailed structure of this part of the star has been the 
subject of a considerable theoretical effort (. 1 > 2 ' 3 ' 4 ^ 5 > 6 > 7 and references therein). At 
the bottom of the crust the appearance of non-spherical nuclei, forming exotic struc- 
tures, has been predicted. In this region besides the Coulomb interaction also the 
quantum effects, associated with a neutron scattering on nuclear inhomogeneities, 
play an important role, leading eventually to a disordered phase . 8 ' 6 . 

In order to understand various processes associated with thermal evolution of 
a neutron star, the properties of Coulomb crystals in the inner crust have been 
studied by many authors (see e.g. . 9 and references therein). The plasma frequency 
of the system reads: uj p — y/ Aitpi on Z 2 e 2 /M , where pi on denotes the ion density, Z 
is the proton number of an ion, and M is the ion mass. Since pi on — 1/ |-7ri? 3 , where 
R c is the Wigner-Seitz cell radius then: 



The quantities: Z and R Cl in the above equation, can be obtained from calcula- 
tions based on a density functional method. However, the determination of the ion 
mass M requires a certain caution, since the nucleus is immersed in the fermionic 
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environment and its bare mass has to be renormalizcd to take into account the in- 
teraction with surrounding neutrons. If the neutrons outside a nucleus have formed 
a normal Fermi gas the problem would have been rather complicated. 10 . In the inner 
crust, however, it simplifies due to the fact that at low temperatures, T < 0.1 MeV, 
the neutron gas is predicted to be strongly paired. In this case one may apply the 
hydrodynamic approximation which describe the system by an irrotational veloc- 
ity field. Although at finite temperatures there always exists a gas of quasiparticle 
excitations which form a normal flow, at temperatures much lower than the critical 
temperature its influence can be neglected. 

We assume in the following that the nucleus has a uniform nucleon density pi n 
and fills a sphere of radius R. The neutron gas is characterized by the density p out . 
We assume moreover that the flow velocity is below the critical velocity for loss of 
superfluidity. 

In order to calculate the mass M of an ion, the Poisson equation for the velocity 
field has to be solved, both inside and outside a nucleus. The boundary conditions 
read: 

&in\r=R = &out\r=R, (2) 

d d 

Pin(-^^in ~ ft ■ u)\ r = R = /w( — $ OM t - H ■ u)\ r =R, (3) 
$out|r-KX> = 0, (4) 

where and <fr out stand for the velocity field inside and outside a nucleus, respec- 
tively. The vector n is an outward normal to the surface of a spherical nucleus and 
r is a radial coordinate (the center of a nucleus is placed at r = 0). The velocity of 
the surface element was denoted by u. The first equation stems from the require- 
ment that the phase of the wave function of the supcrfluid system is continuous, 
the second one is a conservation of mass in radial flow, and finally the third one 
ensures that the correct asymptotic behavior for the velocity field is obtained. The 
knowledge of the velocity field allows us to determine the kinetic energy and corre- 
sponding effective ion mass. Consequently, the following expression for the mass of 
an ion being a subject of a translational motion can be derived: 

4 (1 - -y) 2 

M r ™ = - n & mpin i—lL, (5) 

where 7 = Pout/Pin, and m is a nucleon mass. 

The above expression relates the mass of an ion to the density of the neutron 
gas. Note, that for a uniform system, i.e. 7=1 the ion mass is equal to zero. 
Since the nucleon density vary substantially across the inner crust, one expects 
that the properties of the Coulomb crystal will be noticeably modified. Indeed, the 
plasma frequency in the inner crust changes as compared to the value calculated 

for the case: p ou t — (ion in the vacuum). Namely, — = \ ^ where wl en 

corresponds to the plasma frequency for the system with renormalizcd ion masses. 
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Fig. 1. Plasma frequency in the inner crust as a function of the nucleon density. Nuclear parameters 
were taken from , 4 ' 11 . 




Fig. 2. The reduced lattice thermal energy: u = U/NkgT (left figure) and specific heat per unit 
volume (right figure) at the temperature T = 0.1 MeV, as a function of the nucleon density for a 
bec Coulomb crystal. N is the number of ions and kg denotes the Boltzmann constant. 



In the Fig. 1 we have plotted fiw v and hjS^ n as a function of the nucleon density. 
One can notice large discrepancies between the value obtained using the bare nu- 
clear mass and the renormalized one. It clearly indicates that the Coulomb crystal 
becomes more stiff at the bottom of the crust a . Consequently the density of the 
phonon spectrum is also modified as compared to the case of bare ion masses. The 
thermodynamic functions: the lattice thermal energy and the specific heat have 
been shown in the Fig. 2. One can notice that the renormalized values at a given 

a One expects the phase transition to non-spherical nuclear phases to occur at the bottom of the 
crust. We disregard this effect in this paper. 
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density and temperature are smaller. In the limit of T — > the analytic expres- 
sions for the ratios of lattice thermal energy U and specific heat C v can be derived. 
Applying the phenomenological relations from the ref. 9 one gets: 



where U ren and C£ e " denote the renormalized values of lattice energy and specific 
heat, respectively. 

The following conclusions can be drawn from the presented results: 

(1) There is a substantial renormalization effect of a nuclear/ion mass in the inner 
crust of a neutron star, due to the presence of a superfluid neutron liquid. 

(2) Consequently, the phonon spectrum and thermodynamic functions of the 
Coulomb crystal are significantly altered. 

(3) Thermal and electric conductivities of the inner crust, governed by the electron- 
phonon scattering are expected to be modified. In particular, the contributions 
coming from Umklapp processes have to be recalculated using the renormalized 
ion masses. 

(4) There is clearly a need for consideration of other nuclear degrees of freedom 
like e.g. shape vibrations, which may significantly modify the electron transport 
properties. 
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